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Abstract. We summarize a theoretical framework based on global time- 
reparametrization invariance that explains the origin of dynamic fluctuations in glassy 
systems. We introduce the main ideas without getting into much technical details. 
We describe a number of consequences arising from this scenario that can be tested 
numerically and experimentally distinguishing those that can also be explained by other 
mechanisms from the ones that we believe, are special to our proposal. We support our 
claims by presenting some numerical checks performed on the 3d Edwards-Anderson 
spin-glass. Finally, we discuss up to which extent these ideas apply to super-cooled 
liquids that have been studied in much more detail up to present. 
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1. Why glasses? vs. universality in glassy dynamics 

It is common to encounter in nature systems that resist equilibration with their 
environments and display what is called glassiness. The name is derived from what 
we normally know as glasses, an irregular array of silicon and oxygen atoms without 
crystalline order, much as a liquid, but as hard as a solid. The molecular diffusion 
within glasses is extremely sluggish, slowing down by over 10 orders of magnitude as 
the temperature is slightly decreased near the operationally defined glass transition 
temperature. Hence, the term glassiness became generically associated with very slow 
dynamics p]. 

So, why glasses? The answer to this question has been the focus of much research 
effort for long. It is certainly a rather difficult question, and there have been a number 
of ideas lined up for trying to understand how material systems become glassy. It is not 
even clear whether in many systems there is a thermodynamic phase one can call glass, 
or there is simply a dynamical crossover at low enough temperatures. 

Do we need to fully answer why glasses? before we really further our understanding 
of glassy dynamics! We take the point of view that, by starting from the fact that glassy 
systems exist (as nature presents us with concrete examples), we can then attempt to 
characterize whatever possible universal properties there are in glassy dynamics. 
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To make this statement clearer, let us turn to a question that Anderson poses in 
the introduction of his Concepts in solids text [2]: why solids? This question, again, 
is a rather complex one, and it is not untwined from the question why glasses? if we 
focus on why a regular array of atoms, as opposed to a random packing, forms in the 
first place. Even if one assumes that a crystalline structure forms at low temperatures, 
a detailed quantitative analysis of the energetics remains to be done so as to determine 
if the packing is hexagonal, cubic, body or face centered, etc. Nonetheless, if one starts 
from the existence (as observed in nature) of a state with broken translational symmetry, 
one can construct theories of lattice vibrations (and quantize it) and of electronic band 
structure (and discern between insulators, metals and semiconductors). Solid state 
physics starts from the solid. 

The approach we review in this paper relies on a similar philosophy: we do not 
claim any theory of the glass transition, and we do not attempt to answer why glasses? 
with this particular approach. Our theory does not allow us to make non-universal 
predictions, such as what the glass transition temperature is (if one can be defined) for 
a certain material, or whether the material displays glassy behavior at all. We aim at 
understanding if there is a set of principles, guided by symmetry considerations, that 
can be used to understand certain universal aspects of glassy dynamics once the glass 
state is presented. For example, glasses age [3]. We thus expect that there exists a 
unified approach to describe aging phenomena, and we seek some guiding principles, 
based on dynamical symmetries, that could allow us to understand universal properties 
in the aging regime, including the scaling of spatial heterogeneities. 

We propose that the symmetry that captures the universal aging dynamics 
of glassy systems is the invariance of an effective dynamical action under uniform 
reparametrizations of the time scales [l]-[9]. Such type of invariance had been known 
to exist since the early days in which the mean-field equilibrium dynamics of spin- 
glasses was tried to be understood [TOj, [H] and it was later encountered in the better 
formulated out of equilibrium dynamic formalism applied to the same systems [T2l [IB"] . 
The invariance means that in the asymptotic regime of very long times a family of 
solutions to the equations of motion is found. This 'annoyance', we claim, has actually 
a physical meaning and implications in the fluctuating dynamics of real glasses. In 
order to make this statement concrete, it is necessary to show that the invariance 
also exists in finite dimensional glassy models. With this aim we showed that global 
time reparametrization invariance emerges in the long times action of short-range spin- 
glasses assuming causality and a separation of time scales |l]-[6]. Basically, the second 
assumption amounts to start from a glass state, where one may claim there is a 
separation between, roughly, a time regime where relaxation is fast and another where 
relaxation is slow. The invariance can then be used to describe dynamic fluctuations in 
spin-glasses and, we conjectured, in other glassy systems as well. 

Physically, the emergence of time reparametrization invariance can be though of 
in the following way. The out of equilibrium relaxation of glassy systems is well 
characterized by two-time functions, either correlations or linear responses. In the slow 
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and aging regime they depend on two times and time-translation invariance is lost. The 
proper measure of 'time' inside the system is the value of the correlation itself, and not 
the 'wall clock' in the laboratory. For instance, in a spin-glass the proper measure of 
sample age is the spin-spin overlap or in particle systems it is the incoherent correlation 
function. Age measures can fluctuate from point to point in the sample, what we called 
heterogeneous aging, with younger and older pieces (lower and higher values of the 
correlation) coexisting at the same values of the two laboratory times. The fact that 
the effective dynamical action becomes invariant under global time reparametrizations, 
t — > h(t), everywhere in the sample, means that the action weights the fluctuations 
of the proper ages, C(r;ti,t 2 ), directly, and the times t\ and t 2 in the action are just 
integrated over as dummy variables. To draw an analogy, in theories of quantum gravity 
the space-time variables X^r, a) are the proper variables, and the action is invariant 
under conformal transformations of the world-sheet parameters r and a. 

So what does global time-reparametrization invariance symmetry concretely teach 
us about observables in glasses? So far we discussed a global symmetry or invariance 
with respect to uniform time-reparametrizations. By looking at spatially heterogeneous 
reparametrizations, we can predict the behavior of local correlations and linear 
susceptibilities and the relations between them. For example, we predict that, after 
a convenient normalization that we explain in the main text, the triangular relation 
between the local coarse-grained correlations, C(r;ti,t 2 ), C(r;t2,t 3 ) and C(r;ti,t 3 ), as 
a function of the intermediate time t 2 , £3 < t 2 < ti, at all spatial points, r, should be 
identical to the global triangular relation, in the asymptotic limit of very long absolute 
times and delays between them and very large coarse-graining linear length [9]. Different 
sites can be retarded or advanced with respect to the global behaviour but they should all 
have the same overall type of decay. Similarly, the relation between local susceptibilities 
and their associated correlations should be identical all over the sample p2 6J leading to 
a uniform effective temperature |14j . 

The purpose of this article is to describe our current understanding of dynamic 
fluctuations (heterogeneities) in the non-equilibrium relaxation of glassy systems arising 
from the time-reparametrization invariance scenario [1]-[9],[T5]. We illustrate it by 
presenting critical tests. 

The structure of this review is the following. In Sect. [2] we explain the origin of 
the time-reparametrization invariance scenario. We do not present detailed derivations 
that were already published but we aim at highlighting the main ideas behind the scene. 
In Sect. [3] we list several measurable consequences of the theory. We discuss how they 
have been examined numerically in different glassy models. The discuss here which 
consequences could also be explained by other approaches and which, we believe, are 
unique to our scenario. Finally, in Sect. H] we discuss the scenario. Since this research 
project is not closed yet, we present some proposals for numeric and experimental tests 
as well as some ideas about analytic calculations that could help us understanding the 
limitations of our proposal. 
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2. Time reparametrization invariance 

In this Section we explain how global time-reparametrization invariance develops 
asymptotically in the aging regime of glassy models. For the sake of simplicity we focus 
on the classical formalism and at the end of this section we mention the modifications 
introduced by quantum fluctuations. 

2.1. Mean- field models - dynamic equations 

Schematic models of spin-glasses, structural glasses and ferromagnetic clean coarsening 
are encoded in the family of p-spin models defined by [T6] 

H — — ^ ^ J i\ii...i v s i\ s ii ■ ■ ■ s i p (1) 

with quenched disordered exchanges distributed with the Gaussian law P{Ji 1 i 2 ...i p ) oc 
e -P lJ i 1 i 2 ...i p /( 2NP )_ ^he dynamic variables Sj, i = 1, . . . , N, are of Ising type, Sj = ±1, 
or satisfy a global spherical constraint Y2iLi s l = N. p is an integer parameter: with 
p = 2 and Ising variables one mimics spin-glasses, with p > 2 and Ising or spherical 
spins the phenomenology of structural glasses is recovered, and with p = 2 and spherical 
spins one describes ferromagnetic domain-growth in clean systems. The sum runs over 
all p-uplets; for this reason these models are 'mean-field' in the sense that the saddle- 
point evaluation of the partition function or the dynamic generating functional is exact 
in the thermodynamic limit, N — > oo. The Hamiltonian ([1]) also represents the potential 
energy of a particle with position s — (si, . . . , sn) on an infinite dimensional hypercube 
(si = ±1) or hypersphere (J^i=i s 1 = N). 

Dynamics is introduced with a Langevin equation that represents the coupling of 
the spins to an equilibrated thermal environment. Ising spins are then replaced by 
soft variables by introducing a double- well potential energy, Y^h=i V( s i)i with V(sj) = 
a(s 2 — l) 2 . The initial condition is usually chosen to be random thus mimicking a rapid 
quench from infinite temperature to the working temperature T. In the limit iV — > oo 
exact Schwinger-Dyson equations couple the global correlation and instantaneous linear 
response 

N N 



NC{t,t w ) = J2^{t)si{t w ) , NR(t,t w ) = J2 ' )S,(I) 



(2) 

h=0 



5hi(t 

(X^ili s i(t) = f° r a ^ t)' The field hi couples linearly to the spin variables and, in 
general, we are interested in perturbing fields that are uncorrelated with the equilibrium 
states of the systems. It is not necessary to average over thermal noise or quenched 
disorder since these quantities do not fluctuate in the out of equilibrium regime reached 
when the infinite volume limit (N — > oo) has been taken at the outset. Here are in 
what follows times are measured from an origin that corresponds to the quench to the 
working temperature. Our study applies then in the order of limits 

lim lim (3) 

t w — >oo N— >oo 
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in which the exact causal Schwinger-Dyson equations for spherical models at times t >t w 
read [El US] 

ft rt m 
(d t - z t )C(t, t w ) = / dt' £(t, t')C(t', t w ) + / dt' D(t, t')R(t w , , (4) 

(d t -z t )R(t,t w )= I dt' E(t,t')R(t',t w ) + S(t-t w ) , (5) 
where the vertex, D, and self-energy, S, are functions of C and R 

D(t,t w ) = ^C p -\t,t w ) , E{t,t w ) = P^plc p - 2 (t,t w )R(t,t w ) . (6) 

The Lagrange multiplier z t is fixed by requiring C(t,t) = 1. For Ising problems the soft 
spin constraint can be treated in the mode-coupling approximation [17] or else, one can 
apply a T — T c expansion taking advantage of the fact that the phase transition is of 
second order for p = 2 (Sherrington-Kirkpatrick or SK model) [T3j . 

In such fully-connected models all higher order correlations and linear responses 
factorize and can be written as functions of these two-time functions. Self-consistent 
approximate treatments of interacting particle models with realistic potentials, like the 
mode-coupling approach, yield similar equations with the addition of a wave-vector 
dependence that in the present context can also be taken into account by considering 
models of ci-dimensional random manifolds embedded in N —>■ oo dimensional 
spaces HB1[T2]. 

2.2. Structural glasses: the p > 3 cases 

We now focus on the p > 3 cases that mimic the structural glass problem. We mention 
at the end of this subsection the Ising p = 2 (SK) spin-glass case that is not conceptually 
different but just technically more involved. In Sect. 12.51 we discuss the p = 2 spherical 
problem that yields a mean-field description of coarsening phenomena and is rather 
different from the point of view of time-reparametrization transformations. 

Equations (HJ and ([Hare causal and be solved numerically by constructing C(t,t w ) 
and R(t,t w ) from the initial instant t — t w — 0. An analytic solution is possible in the 
asymptotic limit, as we discuss below. We first present the main features of C and R 
and we later explain how these are obtained from the asymptotic analytic solution. 

Equations (jlj) and ([5]) have a dynamic transition at a critical temperature T^p). At 
T > Td the dynamics occurs in equilibrium and close to the decay of the correlations 
slows down as in super-cooled liquids with the a relaxation time diverging as a power law 
of T — T d [161 120] ■ Below T d eqs. ^ and ([5]) admit a unique solution [121 [13] that is no 
longer stationary. The behaviour of the low-temperature correlation and susceptibility 
is sketched in Fig. [TJ 

The low-temperature solution presents a separation of time-scales. In the long 
waiting-time limit the self-correlation and integrated linear response or susceptibility, 
x(t,t w ) = f* dt' R(t,t'), can be written as 

C(t,t w ) = C s t(t — t w ) + C a g(t,t w ) , (7) 



CONTENTS 



o le-01 




1e-02 



1e-01 



1e+05 



1e+00 



1e-01 



aging & slow (% aa ) 




rapid & stationary (% st ) 



1e-01 



1e+01 



1e+03 
t-tw 



1e+05 



Figure 1. Sketch of the relaxation of the self-correlation and susceptibility in the 
glassy regime. The separation of time-scales is clear in the figure. The Edwards- 
Anderson parameter, q ea the corresponding susceptibility Xea and the asymptotic value 
limt—xxj \(t, t w ) are indicated with horizontal lines. 



X(t, tw) = Xstit - t w ) + Xag(t, t w ) . (8) 

The first terms in the right- hand- side describe the stationary regime at short time- 
differences in which the correlation and susceptibility relatively rapidly approach a 
plateau at limt-^-oo lim^oo C(t, t w ) = q ea and limt-^^co lim^— x(t, t w ) = (1 - 
Qea)/T = Xea- The second terms are the aging relaxation of C towards zero (in 
the absence of an external field), and the aging response of the system towards the 
asymptotic value Xea + Qea/T e ff with T e ff a parameter with the interpretation of an 
effective temperature [H]. 

The stationary and aging relaxation are fast and slow in the sense that 

d t C st (t, t w ) ~ C st {t, t w ) C>q ea , (9) 

d t C ag (t, t w ) < C ag {t, t w ) C <q ea , (10) 

and 

d t Xst(t,t w ) ~ Xst(t,t w ) X<Xea, (11) 

9tXag(t,t w ) < Xag(t,t w ) X > Xea ■ (12) 

The aging self-correlation and susceptibility scale as 

C ag {t, t w ) « q ea f c (^y) , Xa 9 {t, t w ) w q ea f x (^yj • (13) 

The scaling functions satisfy the limit conditions / c (l) = 1, f c (°°) — 0, / x (l) = 
and f x (oo) = 1/T e //. Using mathematical properties of monotonic two-time functions 
one can show that such a scaling holds asymptotically in each (two) time-scale of the 
evolution [13]. While in a system undergoing finite dimensional coarsening TZ(t) has 
a natural interpretation as the typical domain radius, in mean-field models there is 
no immediate understanding of the 'clock' TZ(t) that, in a sense, sets the macroscopic 
time-scale. The numerical solution suggests that TZ(t) is just a power of time. 
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In the asymptotic limit in which the additive separation of time-scales with the 
scaling form ffl3|) holds it is convenient to use a parametric description of the dynamics 
in which times do not appear explicitly. More precisely, the approach to the asymptotic 
scaling, and f c and f x , can be put to the test by constructing 'triangular relations' 
between correlations and susceptibilities, respectively. For generic three long times 
ti > t 2 > t 3 ^ to one computes the correlations (7(2^, 11 > v — 1,2,3. If the times 
are such that the ratios TZ(t^)/TZ(t u ) remain finite in the asymptotic limit, that is to 
say C{t^t v ) = Cag{t^,t u ), one has 

C(*l,t 3 ) = Qea fc {^[Cih^/qea] J ^ \C % , t 3 ) / ' q ea ]} . (14) 

If, instead, t\ = t 2 + t with r > finite, and lZ(t 2 )/lZ (£3) finite in such a way that 
C(ti,h) = (lea + CstiU - t 2 ) and C(t 2 ,t 3 ) = C ag (t 2 ,t 3 ) one has 

C(t 1 ,t 3 ) = mm[C(t 1 ,t 2 ),C(h,t s )] (15) 

asymptotically. This form goes under the name of dynamic ultrametricity. In the 
opposite case t 3 = t 2 — r and lZ(ti)/lZ(t 2 ) finite dynamic ultrametricity also holds. 
These relations follow immediately from the additive separation of time-scales ([8]) and 
the scaling (fl3l) but they can be shown without assuming dynamic scaling, just by 
using the monotonicity properties of temporal correlations [1 3j . The simplest way to 
see these relations at work is to display C(ti,t 2 ) against C(t 2 ,t 3 ), for a chosen value of 
C(ti,t 3 ) < q ea , in a parametric plot in which t 2 varies from t 3 to ti. In the asymptotic 
limit t 3 — > 00 the construction reaches a stable master curve as displayed in Fig. [2]- 
left. The vertical and horizontal parts correspond to t 2 such that C(ti,t 2 ) > q ea and 
C(t 2 ,t 3 ) > q ea , respectively, and dynamic ultrametricity holds. The curved part is for 
t 2 such that all correlations are in the aging regime and its functional form is fully 
determined by f c . The 'clock' 1Z yields the speed at which the data-point moves on the 
parametric curve. A similar construction can be done for the susceptibilities. 

The stationary correlation, C st , and susceptibility, Xst, are linked by the equilibrium 
fluctuation dissipation theorem (fdt), \st = (1 — C st )/T. In the aging regime, instead, 
there is a non-trivial relation between x an d C: Xag = (Q ea — C ag )/T e ff that yields 
f x {x) = (l—f c (x))/T e ff. This relation is also better appreciated if shown in a parametric 
construction in which times do not appear explicitly. In the long waiting-time limit the 
plot x(t, t w ) against C(t, t w ) with t the parameter running from t w to infinity approaches 
a broken line form with the slopes — 1/T (for C > q ea ) and —1/T e ff (for C < q ea ). Again, 
the 'clock' 71 yields the speed at which this curve is constructed upon increasing t. 

An analytic solution to the Schwinger- Dyson equations was derived in the limit 
of long waiting-time in which the separation of time-scales, that is to say the plateaus 
in C and \i are fully established. In the aging regime, one uses the fact that the 
variation of the correlation and linear susceptibility are negligible with respect to all 
terms in the right-hand-side of the equations and can thus be dropped. Furthermore, 
one approximates the integrals by separating the contributions from the stationary and 
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Figure 2. Sketch of the parametric representation of the correlation and susceptibility. 
Left: triangular relation between the correlation function in the asymptotic limit 
£3 — > 00. The three curves correspond to different tx'a such that C(ti,t 3 ) takes three 
values. The breaking points lie at q ea . The arrow indicates the sense of the evolution 
when t% increases from t% to t\. Right: susceptibility, X^i^'i) against correlation, 
Gitii ts) at fixed £3 using > t% as a parameter in the long t% limit. The breaking point 
at (<7ea, Xea) separates the stationary regime where the equilibrium fdt is satisfied from 
the aging regime where it is modified. The slopes are — 1/T and — 1/T e //, respectively. 
The arrow also indicates here the sense of the evolution when t<i increases from £3. 

aging regimes [15] . Equation (jSJ) becomes 

HocRagit, t w ) ~ P(P ~ ^ f dt' Cl g \t, t')R ag (t, f) R ag (t', t w ) (16) 

(fioo is a constant with contributions from lim^oo z t and border terms in the integrals). 
The companion eq. (jl]) takes a similar form within the same approximation. Now, the 
surprise is that the approximate equations are invariant under the transformation 

t — > h t = hit) \ ^ a 9^^ w ' ~^ C ag (h t , h tw ) , 

[ Ra 9 {t,t w ) — > h tw R ag (h t , h tw ) , 

with /if positive and monotonic and = dh tw /dt w . While the functions f c and 
f x and their FD relation are fixed by the remaining approximate equation (fT6|) and 
its companion, the time-reparametrization invariance does not allow one to compute, 
analytically, the clock TZ(t). This problem is similar to the velocity selection problem 
present for instance in Fisher differential equation describing front propagation and the 
like [21J. The exact Schwinger- Dyson equations do have a unique solution with a special 
function TZ(t) that is selected by the short-time difference effect of the time-derivatives. 
However, as time increases and time-differences increase too the effect of the time- 
derivatives diminishes. In the approximate analytic solution we take advantage of this 
fact to solve the equations asymptotically but we introduce in this way a symmetry that 
does not allow us to fix 1Z(t). We obtain, instead, a family of solutions parametrized by 
h t . It is important to reckon that the parametric constructions in Fig. [2] are independent 
of the clock and thus are fully determined by the approximate treatment. 



CONTENTS 



10 



The case p = 2 with Ising spins or Sherrington-Kirkpatrick model has a more 
complicated scaling form with a sequence of two-time scales leading to dynamic 
ultrametricity for all C < q ea asymptotically [13]. This behaviour is technically more 
involved but, as far as the symmetry properties are concerned, it is similar to the case 
treated above. The full dynamic equations have a unique solution but the approximate 
ones acquire time-reparametrization invariance. We shall not discuss these cases further 
in the rest of this review. 

2.3. Short-ranged models - dynamic action 

In a series of papers [l]-[9] we claimed that the time-reparametrization invariance thus 
far introduced via the asymptotic solution to the dynamic equations in mean-field glassy 
problems is indeed an asymptotic property of the dynamics of glassy systems, mean-field 
and finite dimensional. The separation of time-scales stationary-aging has been observed 
in a variety of glassy systems with numerical simulations and experiments [22]. The 
slowness of the decay in the aging regime, eqs. (fTUj) and (fl2"l) and a weak long-term 
memory of the kind (113p . are the hallmark of glassy relaxation. The idea is then that 
time-reparametrization invariance is the symmetry associated to the dominant dynamic 
fluctuations in these sytems. 

In order to pursue this idea forward one has to first prove that the symmetry of 
the saddle-point equations is also a symmetry of the action in the dynamic generating 
functional not only of fully-connected spin models of the mean-field type but also of finite 
dimensional glassy systems. In [I] we derived and studied the symmetry properties of 
the dynamic action - the so-called Martin-Siggia-Rose action associated to Langevin 
stochastic dynamics - of the disorder averaged soft-spin 3d Edwards-Anderson (ea) 
model of spin-glasses. (H = Yluj) Jij s i s j with Gaussian random variables with zero 
mean and taking non-zero value only on nearest neighbours on the (i-dimensional lattice 
and Si = ±1.) In our analysis we took a number of steps that we briefly recall here. 
First, we introduced four fluctuating two-time fields defined on the lattice sites, 

Qf{t,t w ) , with i = 1,...,N and a, 6 = 0, 1 . (18) 

Their thermal averages are the expected values of the local two-time self-correlation 
(a = b = 0), the retarded linear response (a = 0, b = 1), the advanced linear response 
(a = 1, b = 0), and a fourth observable (a = b = 1) that vanishes if causality is preserved. 
Second, we assumed that a separation of time-scales fast-slow, of the type described in 
the previous subsection, applies to these fluctuating fields too. Third, we determined the 
long-time action by using a Renormalization Group (RG) scheme in the time variables. 
This allowed us to write the full action as a sum of two contributions: one from the fast 
regime holding at short time-differences, another one from the slow regime valid at long 
time-differences. The coupling between these two vanishes asymptotically. Fourth, we 
analyzed the surviving terms in the action, that are just the slow contribution. Using 
advanced and retarded scaling dimensions that are just the labels a and b of the fields, 
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one finds that the global time-reparametrization, 

t^h = hit) , Qf(t,t w ) = (h t r(h tw ) b Qf(h t ,h tw ) , (19) 

leaves all surviving terms unmodified. This step is concisely carried out as follows. Take 
a generic term in the action. Under (|T9~I) it transforms as 

Y[dt v ► / ]]dt u (h tv ) A » ■■■ 

J u=l J v=l 

where M is the number of time integrals, the dots represent a product of the fields 
Qf 3 and A u is the sum of all factors h tv arising from the transformation of the fields. 
Interestingly enough, one finds that all the A u equal one. Thus, with a simple change 
of variables one absorbs each factor h tv in the corresponding integration variable and 




Note that in order to prove invariance under the simpler and more common scale 
transformation, h t = Xt, it is enough to have Y^f=i ^ = N. Scale invariance is included 
in the larger global time-reparametrization symmetry but it is, clearly, more restrictive. 
Finally, we showed that the measure in the functional integral is also global time- 
reparametrization invariant, completing the program. We refer the reader to Refs. p[] 
and [6J for the technical details leading to these results. 

In the disordered 3d EA model studied we carried out the disorder average. The 
presence of quenched disorder gave us an analytic control of the theory but this does 
not necessarily mean that such a symmetry develops only for the long time regime 
of models with quenched disorder. It appears that if the essential assumptions are 
causality and unitarity, and a separation of time scales that takes the action to a non- 
trivial asymptotic state (some glassy state), then one expects the symmetry to exist 
for systems without quenched disorder but that are glassy nevertheless. In order to 
check the development of time-reparametrization invariance in problems of particles in 
interaction in finite dimensions one should first obtain the relevant action to work with. 
A good candidate for a starting point is the Dean-Kawasaki stochastic equation for the 
evolution of the local density [23J. The idea is then to write its dynamic generating 
functional and study the symmetry properties of the effective action assuming that a 
separation of time-scales exists. We are currently carrying out this study. 



2.4- Turning a nuisance into something useful - symmetry as a guideline 

The global time-reparametrization invariance implies that the action describing the long 
time slow dynamics of a spin-glass is basically a "geometric" random surface theory, with 
the Q's themselves as the natural coordinates. The original two times parametrize the 
surface. Physical quantities, as the bulk integrated response x(ti,h) = dt' R{t' , t 2 ) 
and correlation C(ti,t 2 ) have scaling dimension zero under t — ► h(t) as well as their 
local counterparts. The emergence of this gauge-like symmetry, which may appear first 
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as a nuisance that relates too many solutions for just one problem. However, it may 
provide a simple way to understand spatial fluctuations in systems that possess this 
global time-reparametrization symmetry. 

Here, a simple analogy with the problem of a ferromagnet may elucidate the point 
we want to make. In a ferromagnet the action is invariant under uniform rotations of 
the magnetization vector m. In the ordered phase, rotation symmetry is spontaneously 
broken, and a certain magnetization direction m is picked. Typically, a vanishingly 
small pinning field selects this direction. The low action excitations are the spin waves, 
fluctuations of the uniform, symmetry broken, state. These spin waves can be described 
in terms of smooth spatially fluctuating rotations around the uniform magnetization 
state. The spin waves, generated by using slowly varying local rotations, are the 
Goldstone modes of the ferromagnet. 

Similarly, in the aging regime of the glassy systems, the action has a global 
symmetry, under uniform time-reparametrizations, t — > h(t), with the fields 
transforming as in eq. ( 1191) . The probability weight of having certain local two- 
time correlation and response, the observables, should be independent of this 
reparametrization. After coarse-graining over a linear length I the non-vanishing 
fluctuating fields are the local coarse-grained correlation (a = b = 0) and linear response 



with the sum carried over the spins in the volume V? = £ d centered at f, and a is the 
lattice spacing. The transformation (ITTj) is now restated as 



and it is an asymptotic symmetry of the action for the slow coarse-grained degrees of 
freedom. Indeed, the symmetry breaking terms, that have their origin in the short- 
time dynamics and short-time difference dynamics, are not identical to zero but become 
vanishing small asymptotically. The particular scaling function 1Z(t) selected by the 
system is determined by matching the fast and the slow dynamics. It depends on several 
details - the existence of external forcing, the nature of the microscopic interactions, etc. 
In other words, the fast modes which are absent in the slow dynamics act as symmetry 
breaking fields for the slow modes. 

In analogy with the spin-wave fluctuations in magnetic systems, that are dictated by 
the rotational symmetry, we proposed that the smooth fluctuations in the glassy phase 
can be obtained by studying the slow varying, position dependent reparametrizations 
t — > h(f, t) around the one reparametrization lZ(t) selected by the short-time dynamics. 
In other words, we basically proposed that there are Goldstone modes for the 
glassy action which can be written as slowly varying, spatially inhomogeneous time 
reparametrizations. This suggests that the slow part of the coarse-grained local 



(a = 0, b = 1) 






(21) 
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correlations and susceptibilities should scale as 




Xa 9 {r;t,t w ) ~ /, 




with f c and f x the same functions describing the global correlation and susceptibility, 
respectively [eqs. f|T3|) and ff23l ] and the same function h(r, t) scaling the two-time 
correlation and susceptibility on each site r [H EE E]- The sum rules C ag (t,t w ) = 
V' 1 f d d r C ag (r; t, t w ) and Xag(t, t w ) = V' 1 f d d r Xa 9 (r; t, t w ) apply. 

The reason for this proposal is that the global reparametrization invariance in time 
of the dynamic action in this two-time regime leads to low action excitations (Goldstone 
modes) for smoothly varying spatial fluctuations in the reparametrization of time, but 
not in the external form of the scaling functions. As in a sigma model (for example, to 
describe the ferromagnet), the external functions f c and f x fix the manifold of states, 
and the local time reparametrizations correspond to fluctuations restricted to this fixed 
manifold of states (in the ferromagnet, the tilting of direction but not the magnitude of 
the magnetization vector). 

2. 5. The spherical p = 2 case or mean- field domain growth 

The spherical SK model (p = 2) can be solved exactly by analyzing the Langevin 
equation in the basis of eigenvectors of the random matrix |24j . While the correlation 
has a very similar behaviour to the one of the p > 3 cases (see Fig. [TJ-left) , the 
susceptibility is quite different. One can mention that the stationary and aging regimes 
in the linear response are not so sharply separated in this case. If one uses an additive 
separation as in (jSJ) the aging contribution to the integrated linear response, Xa g (t,t w ), 
vanishes asymptotically. More precisely, 



Importantly enough, even though the inequalities (flOl) and (112j) are still valid, a 
careful inspection of all terms in the Schwinger-Dyson equations shows that they are of 
the same order asymptotically. Moreover, the stationary contribution to the equations 
in the aging regime is not just a constant: the corrections associated to the asymptotic 
approach to the plateau in the correlation cannot be neglected. As a result one cannot 
simply drop the time derivatives and the Schwinger-Dyson equations in the aging regime 
are not time-reparametrization invariant but just scale invariant, that is to say, they are 
unchanged by the transformation t — > h(t) = At, with C and R transforming as in ( TT7T) 
and A a positive constant. A similar mechanism, though even harder to prove, applies 
to the dynamic equations for the correlation and linear response of the non-conserved 
dynamics of the 0(N) ferromagnetic model in the large N limit [8]. 

In line with what we explained above, the effective action for the slow degrees 
of freedom of the p = 2 spherical model and, more generally, the <i-dimensional 
ferromagnetic O(N) model in the large iV limit, are not invariant under global time- 
reparametrizations but only under global rescaling of time, t — ► At. This marks 




(23) 
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an important difference between models with a finite aging response and these quasi 
quadratic models with a vanishing aging response [8]. 

This result is important for a number of reasons: first, it suggests that the 
susceptibility, or even the effective temperature, might be intimately related to the 
symmetry properties of the dynamics and consequently of the fluctuations; second, it 
suggests that the mechanism for fluctuations in coarsening systems might be different 
from the one of other glassy problems with finite and well-defined effective temperatures. 
In order to justify the latter statement it remains to be checked whether the reduction of 
time-reparametrization invariance to just scale invariance also holds in finite dimensional 
non-field coarsening. 

2.6. Quantum problems 

In the case of quantum models one introduces the effect of dissipation by coupling the 
system to an environment represented, typically, by an infinite ensemble of quantum 
harmonic oscillators. One then uses the Schwinger-Keldysh formalism to write a 
generating functional and from it, in the fully-connected or infinite dimensional cases, 
one derives Schwinger- Dyson equations similar to the ones above. The asymptotic 
analysis of these equations follows the same steps as in the classical limit [25J [T7j (at 
least in the case of a weak coupling to the bath [26]) and the time-reparametrization 
invariance also applies. 

The appearance of an asymptotic invariance under time-reparatrizations in the 
mean-field dynamic equations was related to the reparametrization invariance of the 
replica treatment of the statics of the same models pHl E7j. The latter remains 
rather abstract. Brezin and de Dominicis [27] studied the consequences of twisting the 
reparametrizations in the replica approach. Interestingly enough, this can be simply 
done in a dynamic treatment either by applying shear forces or by applying heat- 
baths with different inherent dynamics to different parts of the system. More precisely, 
using a model with open boundary conditions one could apply a thermal bath with 
a characteristic time-scale on one end and a different thermal bath with a different 
characteristic time-scale on the opposite end and see how a time-reparametrization 
'flow' establishes in the model. 

3. Consequences and tests 

In this Section we discuss how one can put these ideas to the test by presenting a number 
of consequences of global time-reparametrization invariance that are directly measurable 
numerically and experimentally. The properties that we discuss explicitly are: 

1. A growing dynamical correlation length. 

2. Scaling of the PDF of local two-time functions. 

3. Functional form of the PDF of local two-time functions. 
4- Triangular relations between two-time functions. 
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5. Scaling relations for general multi-time functions. 

6. Local fluctuation-dissipation relations. 

7. Infinite susceptibilities. 

In considering these predictions, we shall separate them into two distinct 
classes. The first class contains predictions that are consistent with other theoretical 
scenarios as well as with the presence of reparametrization invariance. Hence, while 
reparametrization invariance leads to these predictions, this class alone cannot be used 
to argue for the role of the symmetry over other mechanisms. The second class, on the 
other hand, contains predictions that are not natural within other frameworks, and to 
the date of this report have no obvious explanation within other frameworks. Properties 
1 and 2 belong to the first class, properties 3-7 to the second. We shall present these 
properties in detail below. 



3.1. Two-time correlation length 

In equilibrium statistical models one defines the static correlation length, £, eq , from the 
spatial decay of the correlation between the fluctuations of the order parameter measured 
at two space points, ([0(f) - ( 0(f) )][</>{?) - {<l>(r'))])\ P _p\ =A ~ A' d+2 ~^ e" A /^, 
where the angular brackets denote an average over the Gibbs-Boltzmann measure. 
£ eg depends on temperature and, in second order phase transitions, it diverges at T c 
leaving only an algebraically decaying correlation. A dynamic equilibrium correlation 
length characterizes the spatial decay of equal-time correlations in the equilibrium 
relaxation of 'usual' systems. Similarly to the static case one defines £(t) via 
([(f)(r,t) - ( 0(f, t) )} [</>(?, t) - (0(r',£))])|-_p| =A ~ e - A /?(*) ; the angular brackets 
indicate here an average over thermal histories and, for simplicity, we omitted the 
algebraic correction to the exponential decay. The average over thermal noise can 
be traded for an integration over the reference space-point r and one then obtains 
i{t) from V~ l jd d r 6<i>(f,t)6<j>(?,t) \ ]r _p ]=A ~ e~ A ^ where 6</>(r,t) = <j>(r,t) - 
V^ 1 J d d r" (f)(r",t). The correlation length, £(£), depends now on temperature and 
total time. 

In systems with slow dynamics in which the order parameter is a two-time entity, 
a two-time correlation length can be defined in analogy to what we described in the 
previous paragraph: 

W, t)<f>(f, t w ) - ( 0(r, t)<f>(f, t w ) )] 

x [(/>(?, t)<f>(P, t w ) - ( <t>(P, t)<f>(P, t w ) )] \ ~ e -A/«t,t») . (24 ) 

/ \r-r'\=A 

Once again by trading the thermal average by a spatial average ( 0(r, t)<f>(r, t w ) ) becomes 
the global correlation C(t,t w ) and £(t,t w ) is derived from the four-point correlation 

C 4 (A; t, t w ) = i / d d r 5W, t)0(f, t w )]5[<j>(?, t)cf)(P, t w )\ , (25) 

V J \f-r'\=A 
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Figure 3. The correlation length in the M EA model at T = 0.6 < T c and L = 100. 
(a) As a function of t — t w ; (b) as a function of 1 — C; (c) in the scaling form £^ a £ 
against 1 — C. These results are taken from [Hj . The vertical arrow in panels (b) and 
(c) indicates the value of q ea - 



with 5[4>(r, t)(f>(r, t w )] = (p(r,t)(p(r,t w ) — C(t,t w ). The quantity C^(A;t,t w ) measures 
the probability that a fluctuation of the two-time composite field 4>(f,t)<p(f,t w ) with 
respect to its global average C(t, t w ) in the spatial position r affects a fluctuation of the 
same composite field at a different site r' located at a distance A from r and averaged 
over the whole ensemble of reference points r in the sample. 

The numerical analysis of C4 in the low temperature out of equilibrium dynamics 
of the 3d EA model [6j [9] (see Fig. [3]), soft sphere [28] and Lennard- Jones [29] mixtures 
yield 

£(£ t ) ~ < ^ st ^ ~ ^ ^ > ^ ea with 

I £,ag(t,tw) C < q ea 

£ag{t,t w ) ~ #(C) an d a a small power (26) 

(a logarithmic growth is also possible). g(C) is a monotonically decreasing function of C. 
This is a two-time monotonically growing function even at time-lags that are longer than 
the waiting-time dependent a relaxation time. Note the difference with what is observed 
in the super-cooled phase where a number of numerical and experimental measurements 
point at a dynamic correlation length that reaches a maximum at the a relaxation time 
and later diminishes to zero [30]. The divergence has a clear interpretation within the 
global time-reparametrization scenario: it is due to the generation and development of 
the zero mode. The global reparametrization invariance symmetry develops only in the 
limit of very long times, so that at intermediate times the irrelevant terms that are scale 
down to zero are still manifest. These irrelevant terms, symmetry breaking ones, give 
a finite length scale (or a finite 'mass') to the soft reparametrization modes. Because 
they are irrelevant, the correlation length increases or, equivalently, the mass decreases 
asymptotically. 

Even though the numerically accessible times are sufficiently long so as to see 
the separation of time-scales in the relaxation of the global correlation, the correlation 
length is still very short. In numerical simulations £ reaches of the order of 4 lattice 
spacings or inter-particle distances in the spin-glass problem or soft sphere and Lennard- 
Jones system, respectively. £ just increases by, say, a factor 4 when the waiting-time is 
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increased by nearly 4 orders of magnitude. The expected limit £(t,t w ) — > oo is thus far 
from being attained. 

As we already mentioned in the introduction to this Section, the growing length 
scale is consistent with the development of time-reparametrization invariance, but it 
is also consistent with other mechanisms. Hence, it alone cannot be used to argue 
unambiguously in favor of the symmetry-based approach. For example, theories based 
on the mode-coupling approach or random first order scenario [33J and its refinement 
including the effect of entropic droplets |34j . dynamical criticality controlled by a zero- 
temperature critical point [35J and frustration limited domains [36J are used to justify 
the growth of a dynamic length-scale, at least in the super-cooled liquid. Hence, the 
existence of a growing length scale belongs to the first class of predictions we mentioned 
in the introduction to this section. 

Finally, let us note that the fluctuations in the susceptibility, or in multi-time 
correlations, see eq. (|35|) . and associated susceptibilities, can be used to derive other 
correlation lengths. It would be interesting to check whether all these behave in the 
same manner. 

3.2. Scaling of the PDF of local two-time functions 

The most direct way of testing the mere existence of local fluctuations is to measure the 
probability distribution function (pdf) of local coarse-grained correlators, C(f;t,t w ), 
and linear susceptibilities, x(r;t, t„,), at different pairs of times, t and t w . In such 
a measurement one is forced to use finite coarse-graining lengths, i then becomes a 
parameters that has to be taken into account in the scaling analysis of the results. 

In Ref. [9] we showed that, quite generally, the PDF of local coarse-grained 
correlators can be scaled onto universal curves as long as the global correlation, C(t, t w ), 
is the same, and the ratio of the coarse graining length over the dynamical correlation 
length, £/£(t,t w ), is held fixed (see [31] for a similar discussion applied to the super- 
cooled liquid). Such scaling can be easily understood as follows. At fixed temperature 
the PDF p[C r ; t, t w , £, L] depends on four parameters: two times, t and t w , and two 
lengths, the coarse-graining length, £, and the size of the system, L. As in the aging 
regime C(t,t w ) is a monotonic function of the two times and £(t,t w ) ~ t^g(C), one 
can trade the two times by C and £ in complete generality. The next step is a scaling 
assumption: that in the long times limit the PDFs depend on the coarse graining length 
£, the total size L and the scale £ only through the ratios £/£ and £/L. This last 
step, we should stress, is really a scaling assumption, and not a trivial requirement 
from dimensional analysis. The lengths £, L and £ are already dimensionless as they 
are measured in units of the lattice spacing. The end result from the rewriting of the 
parameters in terms of the global correlation and the scaling hypothesis is that the 
PDFs characterizing the heterogeneous constant temperature aging of the system can be 
written as 




(27) 
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Figure 4. PDF of local coarse-grained correlations C r at different times t and t w in 
the 3d Edwards- Anderson model with L = 100 at T = 0.6 < T c . The waiting-times 
are given in the key and the global correlation is fixed to C = 0.4 < q ea . (a) The 
coarse-graining boxes have linear size I = 9 in all cases. The curves do not collapse, a 
slow drift with increasing t w is clear in the figure, (b) Variable coarse-graining length 
I chosen so as to held £/£ approximately constant. The collapse improves considerably 
with respect to panel (a). These results are taken from [9]. 

In Fig. H] we show the scaling of the distribution of local coarse-grained correlations 
in the 3c? EA model and the effect of the scaling variable £/£ (the size of the system, L, 
is sufficiently large so that £/L vanishes in practice). It is noteworthy that a reasonable 
scaling with the global correlation held fixed and not taking into account the effect of 
the second scaling variable has been already achieved, approximately, in the Edwards- 
Anderson model [3 E], as well as in the kinetically constrained models studied in [7] 
and Lennard- Jones systems [32]. This is justified by the fact that £ varies very slowly 
with t w . However, it is clear that at long though finite times one has to hold the ratio 
£/£ constant to obtain a full collapse in all cases. 

The scaling variable £/£ allows one to study the change in the functional form of the 
PDFs upon modifying the coarse-graining volume. Indeed, the pdfs should crossover 
from a non-trivial form to a simple Gaussian when £ goes through the value £. In 
summary, for finite £ one identifies three ^-dependent regimes with different functional 
forms of the PDFs: 

• £ <C £. For finite £ this means £ of the order of the lattice spacing, £ ~ a. In this 
case the PDFs do not have any particular structure. 

• I <~ £. For finite but large £ this case is non-trivial and indeed the one that is 
accessed with numerical simulations and experiments. One finds that the statistics 
is non-Gaussian for all C. The skewness decreases from zero at C — 1 to 
reach a minimum and then increase again at small values of C. As regards the 
functional form, one finds that a Gumbel-like functional form, characterized by a 
real parameter that depends on C and £/£ describes the data rather well for, say, 
q ea /2 ~ C ~ q ea - 

• £ ^> £. In this limit one matches the central-limit theorem conditions and the 
statistics becomes Gaussian for all C. 
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The analysis of the ^-dependent pdfs thus provides an independent way to estimate £. 
3.3. Effective action for local ages 

The argument leading to eq. (127}) is a scaling hypothesis and it does not rely on time- 
reparametrization invariance. Indeed, there exist models in which the scaling (128]) for 
the PDF of local coarse-grained correlations is found, e.g. the O(N) model in the large 
N limit [8], and global time-reparametrization does not apply 

The implications of time-reparametrization invariance appear later, as a prediction 
for the functional form of the asymptotic limit 



In order to study the functional form that can take let us now use the symmetry 
argument to analyse the statistics of the fluctuations of the local correlations. So far 
we have not yet determined how much the h(r, t) vary in space and time. To this end 
we need to derive an effective action for these functions that will tell us how costly it 
is to deviate from the average clock TZ(t). Ideally, one would like to derive this action 
from the microscopic one. This should be possible in quasi mean-field models such as 
the p-spin model with Kac long (but finite) range interactions [37]. For the moment 
we have just proposed the simplest action that serves our purposes in the ideal limit in 
which the zero mode is fully developed £ — > oo and the local quantities are measured in 
the infinite coarse-graining volume limit with £/£ — > [7J. Otherwise the parameter £ 
should be taken into account. 

To start with we worked with the more convenient transformed variable h(r, t) = 
e -<p[r,t) th a t implies 



and we searched for the simplest action that satisfies the constraints due to the 
symmetries. These are: 

i. The action must be invariant under a global time reparametrization t — > h(t). 

ii. If our interest is in short-ranged problems, the action must be written using local 
terms. The action can thus contain products evaluated at a single time and point in 
space of terms such as ip(r, t), d t ip(r, t), V<p(r, t), Vd t ip(r, t), and similar derivatives. 

Hi. The scaling form in eq. (129]) is invariant under <fi(r,t) — > if(f,t) + $(r), with $(f*) 
independent of time. Thus, the action must also have this symmetry. 

iv. The action must be positive definite. 

These requirements largely restrict the possible actions. The one with the smallest 
number of spatial derivatives (most relevant terms) is 



Poo {C r -C) = Km lim lim p[C r - C(t, t w ), £/£(t, t w ), f(t, t w )/L] . (28) 



C(t,t w )=C 




(29) 




K 



d t <f(f,t) 



(30) 
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with K a stiffness. A term M d t tp(f, t) is also allowed by symmetry but since its space- 
time integral is constant we drop it. The action solely depends on the time derivatives 
dtf(r,t) and it is simple to check that it satisfies all the four constraints enumerated 
above (the last requirement follows from the fact that h(f, t) are monotonically 
increasing functions of time) [7J. 

Due to the simple form (130]) the d t tp(r, t) are uncorrelated at any two different times 
ti and t 2 . Thus the expression Atp?^ = dt' d t 'ip(r,t') entering the exponential in 
the scaling form in eq. ( 1291) is a sum of uncorrelated random variables in time. One 
can interpret such expression as the displacement of a random walker with position 
dependent velocities. Alternatively, one can think of the space-dependent differences 
A(ff\\ as the net space-dependent height (labeled by t) of a stack of spatially fluctuating 
layers dt dt<p(r,t). The action for the fluctuating surfaces of each layer is given by 



The statistics of the Ay??^ are completely determined as follows. The action 
( 1301 transforms into one of a Gaussian surface after the introduction of a 'proper' time 
t = \n7Z(t), and the change of variables, ip 2 (r,t) = d T ip(r,r). Indeed, 



Due to the Gaussian statistics of the ijj, it is simple to show that connected iV-point 
correlations of A 
varphi^W^ satisfy 



where the function T can be obtained from Wick's theorem, summing over all graphs 
that visit all sites (connected) with two lines (because of ip 2 ) for each vertex i 
corresponding to a position r«. The reparametrized times appear only in the prefactor 
r(t) — r(t w ) = \\i\Tl{t) /1Z{t w )]. The probabilistic features of the fluctuations of local 
correlations C{f,t,t w ) depend on times only through 7Z(t)/7Z(t w ), and hence only 
through the global correlation itself C(t, t w ). In consequence, the action ([32]) implies 
the scaling (128]) . The fact that the time-dependencies of the statistical properties of the 
two-time local coarse-grained correlations are fully determined by the global correlation 
is a very welcome property of action (|30|) since it was not obvious a priori. 

Having the forms in eqs. (!3T|) and (|32|) allows us to make some quantitative 
predictions about the form of the pdfs. With some algebraic manipulations one shows 
the following generic features [7]: 

• The distribution is non-Gaussian for all C . 

• For C ~ q ea the pdf is negatively skewed and once put into normal form it is 
very close to the distribution of the global equilibrium magnetization in the 2c? XY 
model [38]. It can then be approximately described by a generalized Gumbel form 
with real parameter. 



eq. dHQD. 




(31) 



(32) 



(AwJL ApJ^.-.A^rJ^Jc = [r{t) -r(0] Hn,r 2 , • • • ,r N ) 



(33) 
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• In the opposite limit C ~ the PDF is positively skewed and it does not take any 
recognizable form. 

If one is interested in testing the action further one can simply use eq. (1321) to generate, 
numerically, the ip(r, t), construct the C(f, t, t w ) from these functions, and then compare 
the PDFs thus obtained to the ones measured, say, in a numerical simulation of a given 
problem. Note that the scaling function f c also plays a role in the functional form of 
the PDF of local correlations. The same argument applies to the susceptibilities. 

The local coarse-grained correlations are, by construction, sums of correlated 
random variables (unless I ^> £). With numerical simulations of the 3d EA model [9] 
and kinetically facilitated lattice gases [7j we found that the PDFs of correlations coarse- 
grained over finite lengths i have a functional form that resembles a generalized Gumbel 
distribution characterized by a continuous parameter that depends on and the value 
of the global correlation, C, when C ~ q ea (see also [38]). This fact is consistent with the 
discussion above and also with the observation of Bertin and Clusel that Gumbel-like 
PDFs with real parameter characterize the statistics of sums of random variables with 
particular correlations between the elements [39] . The fact that we obtain Gumbel-like 
distributions then means that the correlations between the terms in the sum are of the 
form needed to get this type of PDF. 

In short, the time-reparametrization scenario predicts, in its simplest setting, that 
eqs. (EH) and ( 1321) fully characterise the fluctuation of the local correlations in the large 
times and coarse-graining volume limits. 

3. 4- Two-time scaling of local functions 

As we argued in Sect. I2.4[ the global time-reparametrization invariance suggests that 
in the ideal asymptotic limit the slow part of the coarse-grained local correlations and 
susceptibilities should scale as in eq. ( 122]) in the ideal limit a <C £ <C £. In practice 
the ideal limit is not reached and one is forced to work with finite correlation lengths 
and thus finite coarse-graining lengths too. The finite I will then play a role and has 
to be taken into account. We now present some tests of eq. (1221) that are based on the 
parametric representation of the dynamics explained in Sect. [2] and take into account 
the finite value of I. 

Let us then imagine that we compute three local coarse-grained two-time 
correlations, CV, at three space points r\, r*2 and r*3, using a given coarse-graining length, 
a <C £, and that we obtain functional forms that are characterized by eq. (122]) with, say, 
h(fi,t) = ln(t/to), h(r2,t) = t/to, and h(fs,t) = e ln2 (*/*°), in the aging regime. In Fig. [5] 
we sketch the decay of these correlations for the same t w as a function of time-delay. 
The plateau is at the same height since q ea as well as the the full stationary decay are 
not expected to fluctuate. The external function f c is the same in all curves. It is clear 
that the decay of the three correlations follows a different pace, the one at r 3 is the 
fastest while the relaxation at f-j is the slowest. 

The simplest way to put the proposal (1221) to the test is to analyze the implications 



CONTENTS 



22 




1e+00 1e+02 1e+04 1e+06 c 
t-tw 

Figure 5. Left: sketch of the decay of the correlation with the same stationary decay 
to q ea - shown with a horizontal dashed line - and three choices of the scaling function 
h(n,t) = ln(i/*o)) h(r 2 ,t) = t/t , and h(r 3 ,t) = e lna (*/*o). The waiting-time is the 
same in all curves. Right: the relation between the integrated linear response against 
the correlation. With a solid line, the parametric plot for fixed and long t w , using t 
as a parameter that increases from t w at C = 1 to oo at C = 0. With symbols, the 
three pairs (Cj{t,t w ),Xj{t,t w )) for the same t w , a fixed value of t and hj(t) as in the 
left panel. 

of eq. (l2"2l on local triangular relations. In Sect. [2] we showed that two-time functions 
with a separation of time-scales as in eq. (|Sj) and an aging scaling as in eq. (fT3|) are 
related in a parametric way in which times disappear, see the sketch in Fig. [21-left. 
Equation (J2"2"j) implies that the local (fluctuating) two-time functions should verify the 
same relation 

C ag (r; h, t 3 ) = q ea f c {f-'lCagir; t u t 2 )/q ea ] f^[C ag {f- 1 2 , t 3 )/q ea ] } . (34) 

This is a result of the fact that time-reparametrization invariance restricts the 
fluctuations to appear only in the local functions h(f,t) while the function f c is locked 
to be the global one everywhere in the sample. 

A pictorial inspection of this relation should take into account the fact that while 
the stationary decay is not expected to fluctuate, the full aging relaxation and, in 
particular, the minimal value of the local two-time functions, C{f]t\,t^), are indeed 
fluctuating quantities. The parametric construction on different spatial regions should 
yield 'parallel translated' curves with respect to the global one, as displayed in Fig. [21- 
left. Fluctuations in the function f c would yield different functional forms in the 
curved part of the parametric construction. A more quantitative analysis can be 
done by using the knowledge of f c that can be extracted from the global correlation 
decay. Indeed, if f c is known, the parametric plot f^ 1 {Cm/ (lea)/ V fc 1 (Cn3/lea) against 
/c~ 1 (Cr23/<Zea)/ a/ f^ 1 (C-m/ Qea) should yield a master curve identical to the global one 
with different sites just being advanced or retarded with respect to the global value. 
This is another way of stating that the sample ages in a heterogeneous manner, with 
some regions being younger (other older) than the global average. (For simplicity 
we used a chort-hand notation, C^ u = C(f;t^t v ) with \i,v = 1,2,3.) If the time- 
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Figure 6. The triangular relation in the 3c? EA model. Upper left panel: the thick 
(black) line represents the global C(t\,t%) against C(t 2 ,t 3 ) using i 2 as a parameter 
varying between t 3 = 5 x 10 4 MCs and t x = 9 X 10 6 MCs, C(t u t 3 ) ~ 0.35 and 
lea ~ 0.8. The curved part is well represented by the hyperbolic form C(ii,ia) ~ 
qeaC(t u t 3 )/C{t 2 ,t 3 ) ~ 0.79 x 0.35/C(* 2 ,*3) that corresponds to / c (x) ~ With 
different points joined with thin lines we show the triangular relations between the 
local coarse-grained correlations on five randomly chosen sites on the lattice (£ = 30). 
Upper right panel and lower left and right panels: 2d projection of the joint probability 
density of C(r; h , *a) > *s)/C(r; h , t 3 ) and C(r;t 2 ,t 3 VC(*i>*3)/C , (r;ti,t3) at 
fixed three values of the intermediate time, t 2 — 1.5 x 10 5 MCs, 8 x 10 5 MCs and 
5 x 10 6 MCs, respectively and I = 10. The global C(ti, t 2 ) against C(t 2 ,t 3 ) using t 2 as 
a parameter is shown with a thick green line. The green points indicate the location of 
C(ti,t 2 ) and C(t 2 ,t 3 ) for the chosen t 2 s. Each point in the scatter plot corresponds 
to a site, r. The lines indicated the boundary surrounding 25%, 50% and 75% of the 
probability density. The cloud extends mostly along the global relation as predicted 
by time-reparametrization invariance. These results are taken from [9]. 
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reparametrization mode is indeed flat the local values should lie all along this master 
curve in the aging regime. 

The conclusions drawn above apply in the strict a <C i <C £ limit. In simulations 
and experiments £ is finite and even rather short. Thus, I is forced to also be a rather 
small parameter, in which case 'finite size' fluctuations in f c are also expected to exist. 
The claim is that the latter should scale down to zero faster (in €) than the fluctuations 
that are related to the zero mode. 

We have tested these claims in the non-equilibrium dynamics of the 3d EA spin- 
glass [9j. The results are shown in Fig. [6j In the upper left panel we show the global 
triangular relation (thick black line) as well as the local one on four chosen sites. The 
separation of time-scales is clear in the plot. The aging part is rather well described by 
f c (x) ~ x~ b and the local curves are quite parallel indeed. In the remaining panels in 
Fig. [6] we show the 2dprojection of the joint probability density of the site fluctuations 
in the local coarse-grained correlations at different chosen times t 2 , t 2 = 1.5 x 10 5 MCs 
(upper right), t 2 = 8 x 10 5 MCs (lower left), and t 2 = 5 x 10 6 MCs (lower right). 
Taking advantage of the fact that f c {x) ~ x~ b we use a very convenient normalization 
in which we multiply the horizontal and vertical axes by [C13/CV13] 1 ^ 2 . Global time- 
reparametrization invariance, expressed in eq. fl34l . implies that the data points should 
spread along the global curve indicated with a thick green line in the figure. Some sites 
could be advanced, others retarded, with respect to the global value - shown with a 
point on the green curve - but all should lie mainly on the same master curve. This is 
indeed quite well reproduced by the simulation data in the three cases, C(ti, t 2 ) close to 
C(ti,t 3 ) (upper right panel), C(ti,t 2 ) close to q ea (lower right panel), and C(ti,t 2 ) far 
from both (lower left panel). Most of the data points tend to follow the master curve 
though some fall away from it. The reason for this is that eq. (1341 should be strictly 
satisfied only in the very large coarse-graining volume limit (£ ^> a) with <C 1 while 
we are here using £ = 10a ~ £, see the discussion in Sect. 13.11 

The triangular relation can be used to test the fluctuations of the local 
susceptibilities too. Indeed, if the separation of time scales (JBJ) and the scaling (fT3l) 
apply to the global susceptibility, the local ones, after the convenient normalization 
by the maximum value, should follow another master curve, identical to the global one. 
Note that time-reparametrization invariance as we discuss it here implies that both local 
correlations and susceptibilities should be fluctuating quantities. 

Finally, notice that, in contrast to the growing correlation length scale, there is 
no blatantly obvious explanation of these triangular relations within other theoretical 
scenarios. These relations are perhaps the most direct consequence of the time 
reparametrization symmetry arguments, and so this prediction falls within the second 
class we discussed in the introduction to this section. 
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3.5. Multi-time scaling 

In general muti-time correlations are non-trivially related to two-time ones. One can 
take as an example a generic coarse-grained connected four-time correlation. If this 
function is monotonic with respect to all times, and the two-point correlations scale as 
in (TT3l) for all pairs (t M , t v ) with /i, v = 1, 2, 3, 4, the four-time correlation should behave 

as 



with the same external function g for all r, in the asymptotic limit in which all times are 
widely separated and the corresponding two-time correlations fall below q ea . Parametric 
constructions could be envisaged to test this relation. 

3.6. Local fluctuation- dissipation relation 

The asymptotic relation between the global correlation and susceptibility 



was first obtained in mean-field disordered models [T2J [13] and later observed in 
simulations of many realistic systems (spins and particles in interaction on finite 
dimensional spaces). — (c?cx(C)) _1 defines an effective temperature [H]. In the aging 
regime, that is to say for C < q ea , three behaviours have been observed in mean-field 
systems: in structural glass models x(C) is linear in C (solid line in the right-panel in 
Fig. [5]); in spin-glass models x(C) is a non-linear function of C; in coarsening systems 
x(C) is a constant equal to (1 — q ea )/T. 

The scaling in eq. (|22|) implies that the parametric construction 'local susceptibility 
against local correlation' should fall on the master curve for the global quantities but 
can be advanced or retarded with respect to the global value; again in the theoretical 
limit a <C i <C This behaviour is sketched in Fig. [fright for the three sites 
whose correlations are displayed in the left panel. The restricted relation between local 
susceptibility and local correlation in eq. (|22|) arise from the fact that the fluctuations 
are due to local reparametrizations alone and not to changes in the external functions 
f c and f x (much as in transverse vs. longitudinal fluctuations in a non-linear cr-model). 

An important property of the interpretation of the fluctuation dissipation relation 
in terms of effective temperatures is that one expects all observables evolving in the same 
time-scale to equilibrate and hence have the same value of the effective temperature [H] 
in an asymptotic regime with slow dynamics and small entropy production. Within 
the time-reparametrization invariance approach the local effective temperature, defined 
from the slope of the x(C) plot, is automatically the same in the whole sample within 
a correlation scale, just because the functions f c and f x do not fluctuate. 

In Fig. [7] we show the joint PDF of local correlations and susceptibilities of the 
3d EA spin-glass in its glassy phase; the accord with the analytic prediction is very 




(35) 




(36) 
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satisfactory [5J E] with the additional spreading away from the master curve ascribed to 
the fact that I is finite and not very different from £. 




Figure 7. (a) The joint PDF p(C r ,Xr) at two times (t w ,t) such that the global 
correlation is C(t, t w ) = 0.7 < q ea in the 3d EA model, (b) Projection of three contour 
levels. The crosses are the parametric construction x(C) for several values of the total 
time t larger than t w . The dotted straight line is fdt at the working temperature T. 
These results are taken from [6]. 



3.7. Infinite susceptibilities 

Zero modes are intimately related to infinite susceptibilities. Indeed, systems with 
continuous symmetries are sensitive to arbitrarily weak perturbations. In the present 
context the approximate global time-reparametrization invariance implies that one 
can easily change the 'clock' lZ(t) characterizing the scaling of the global correlation 
and linear response by applying infinitely weak perturbations that couple to the zero 
mode. An illustration of this property is the fact that the aging relaxation dynamics 
of glassy systems is rendered stationary by a weak perturbing force that does not 
derive from a potential while the x{C) relation in the slow regime is not modified [40J. 
One could envisage more refined tests such as applying a perturbation that imposes 
different scalings on two macroscopic borders of the system and see how the time- 
reparametrization wave develops in the bulk. 

3.8. Conclusions 

In conclusion, the global time-reparametrization invariance scenario gives a mechanism 
for the divergence of the correlation length £ though others have also been proposed in 
the literature. There are a number of predictions, as the parametric relations between 
local coarse-grained correlations measured at different times and the local fluctuation- 
dissipation relations that, to our knowledge, are not explained by other scenarios. As 
regards to the easier to measure pdfs of local correlations the framework is not only 
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consistent with the scaling (1281) - that also arises from simple scaling assumptions - 
but it also justifies the non-Gaussian and Gumbel-like functional form of the pdfs 
that follows from the proposed effective action for local 'ages'. Moreover, systems with 
global time-reparametrization invariance should have as important fluctuations in the 
local susceptibilities as in the local coarse-grained responses. 

4. Discussion 

We presented a summary of studies of dynamical fluctuations in glassy systems that 
are based on the idea that, in the long time regime, a global time-reparametrization 
invariance emerges in the effective action describing the aging dynamics. We discussed 
how this symmetry concretely appears in mean-field spin models, and how it can 
be shown to emerge at the level of the action for short-ranged spin glass models 
with quenched disorder. Two assumptions are used to prove the global time- 
reparametrization of the action for the short-range spin glass model i) causality and 
unitarity, and ii) a separation of time scales between a fast (or stationary) and a slow 
(or aging) relaxation, where in the latter time translation invariance is broken. 

That the dynamical action is symmetric under uniform, i.e., spatially independent 
reparametrizations of time variables (t — > h(t)) suggests that the dynamic fluctuations 
that cost little action should be describable in terms of position dependent, long 
wavelength, reparametrizations of the form t — > h(f,t). These should be the Goldstone 
modes associated with breaking time-reparametrization invariance symmetry. 

We presented predictions of this theoretical framework and tests that we performed 
in the 3d Edwards- Anderson model to falsify these predictions. Among the consequences 
of our theoretical framework are those listed in Sect. [3l Some of them find an 
explanation within other theories as well; others are particularly related to the time 
reparametrization invariance scenario. For example, a correlation length that grows 
in time is associated with the asymptotic approach to the long-time regime in which 
the symmetry is fully developped, and the long wavelength modes eventually become 
massless. The existence of a growing length is also predicted by other models. The 
functional form of the triangular relations relating local two-time correlations between 
three different times is, instead, particular to our framework. 

In this review we showed tests of the predictions of the global time- 
reparametrization invariance scenario performed on a finite dimensional spin-glass 
model [9]. In the past we also studied the dynamics of kinetically facilitated models, 
without quenched disorder, along the same lines [?]. We believe that, if aging dynamics 
is a universal property of glassy systems, then these ideas should also apply to interacting 
particle systems without quenched disorder. The rationale is that the two assumptions 
leading to the global time-reparametrization invariance of the dynamical action, namely 
causality and unitarity, and a separation of time scales, should also hold for other glassy 
systems. The former assumption we can take as a fact. The second is, in a way, the 
assumption that a glassy phase exists, even though we say nothing as of why it does. As 
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we stated in the introduction, we do not aim at the question why glasses?, but instead we 
focus on the possible universal dynamical properties once the glassy state is presented by 
nature. Some of the consequences of the symmetry have already been tested numerically 
in Lennard- Jones systems [291 [32] but there is still much room for more detailed studies, 
including the analysis of local triangular relation between correlation and susceptibilities 
and tests of their joint behaviour. 

We thus propose that the asymptotic global time-reparametrization invariance, 
and the associated low action long wavelength local reparametrizations, constitute the 
mechanism by which dynamic fluctuations, that is to say heterogeneities, are generated 
in the glassy state. Moreover, this mechanism may also apply, in an approximate form, 
to the super-cooled liquid regime. It should just be an approximation because in super- 
cooled liquids the symmetry is not fully developed and there is then a cut-off setting 
the limit of the spatial and temporal extent of the heterogeneities, in sharp constrast 
to the low temperature glassy regime in which the symmetry is realized asymptotically 
and fluctuations of all sizes exist. The growth and divergence of the (two-time) dynamic 
correlation length defined from the study of the space correlation of the two-time order 
parameter is a manifestation of the growth and divergence of these fluctuations in the 
glassy state; in contrast, such growth is interrupted in the super-cooled liquid. 

To better understand the distiction between the glassy state with its asymptotic 
symmetry and the super-cooled state without the fully developing symmetry, consider 
the phenomenology of dynamic fluctuations as a function of temperature. Dynamic 
heterogeneities in the super-cooled liquid phase have been identified numerically and 
experimentally [52J. These are in a number of ways more important than what is 
observed in a simple liquid or a solid. In the super-cooled liquid phase while the 
full relaxation is stationary, there is still a time-scale separation with the correlations 
decaying as a function of time-difference first rapidly to a temperature-independent 
plateau and then slowly towards zero. The latter is the structural or a-relaxation. The 
a relaxation time, t a , is finite in the super-cooled liquid regime and it increases upon 
decreasing temperature. The global parametrization chosen by the symmetry breaking 
terms in the slow regime is TZ(t) oc e - */* a and C(t,t w ) = q e afc(e~^~ tw ^ ta ) in this case. 

The mode-coupling approach to super-cooled liquids is based on approximate 
dynamic equations for the relevant correlators of realistic systems that are very similar 
to the p-spin Schwinger-Dyson equations in the high temperature phase [T6l US] . In 
these equations the correlators are already expressed as functions of the time-difference, 
r = t — t w . Close to the critical temperature the separation of time-scales develops in 
the mode-coupling equations. The approximate analysis of the a relaxation predicted 
by these equations also relies on dropping the r derivatives and approximating the 
integrals by assuming a sharp time-scale separation. The remaining asymptotic (large 
r) equations are invariant under reparametrizations of r. 

We then expect the local coarse-grained correlations and integrated linear responses 
in the super-cooled liquid to be, to a first approximation, stationary (after a sufficiently 
long waiting-time that goes beyond the equilibration time) but with different finite 
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structural relaxation times, fluctuating about the value that characterizes the decay 
of the global correlations. This is consistent with the experimental observation that 
dynamic heterogeneities in supercooled liquids seem to have a lifetime of the order of 
the relaxation time. Deviations from stationarity are not completely excluded for finite 
I but they should be less important than in the aging low-temperature regime. 

There is, however, an important difference with respect to the aging regime, in 
which the equilibration time diverges and local relaxation times or, better stated, local 
ages can fluctuate without limit when t w — > oo. At high temperatures one does 
not expect to find heterogeneities with arbitrary long relaxation time. Furthermore, 
heterogeneities have a finite spatial extent and one can then suppress them by using 
sufficiently large coarse-graining volumes. The correlation length is stationary, £(t), 
and it saturates in the limit of long-time differences, r — > oo. The saturation value, 
though, increases for decreasing temperature. From a theoretical point of view, this 
picture is, in a sense, similar to the one that describes the equilibrium paramagnetic 
phase in the O(N) model, just above the ordering transition temperature. 

When lowering the temperature the size and life-time of the heterogeneities 
increases. A p-spin or mode-coupling approach predicts that their typical size and thus 
linv^oo^r) diverge at the mode coupling transition temperature [33]. In real systems 
the divergence at T c is rendered smooth and £ does not strictly diverge asymptotically. 

At still lower temperatures the bulk quantities age and we expect then to observe 
heterogeneous aging dynamics of the kind described in this review, with a two-time 
dependent correlation length for the local fluctuations. The heterogeneities age as well, 
in a 'dynamic' way. By this we mean that if a region looks older than another one when 
observed on a given time-window, it can reverse its status and look younger than the 
same other region when observed on a different time-window. 

The infinite susceptibility with respect to perturbations that couple to the zero 
mode are illustrated by the fact that the clock of the bulk quantities that is selected 
dynamically is very easy to modify with external perturbations. A small force that does 
not derive from a potential and is applied on every spin in the model renders an aging p 
spin model stationary [ID] while the model maintains a separation of time-scales in which 
the fast scale follows the temperature of the bath, T, while the slow scale is controlled 
by an an effective temperature, T e ff > T. In this case, the aging system selects a 
time-reparametrization 7Z(t) = t while in the perturbed model TZ(t) = e~ t l ta . Similarly, 
the aging of a Lennard- Jones mixture is stopped by an homogeneous shear [41J. A 
different way to modify the time-reparametrization that characterizes the decay of the 
correlations is by using complex thermal baths [42] . 

The picture that we have described applies to long times but not as long as to enter 
the activated regime that we still do not know how to characterize theoretically, not 
even at the bulk level. This regime corresponds to times that scale with the system size. 
The success of mean-field models, or the mode-coupling approach, in describing the 
bulk dynamics of glassy systems, at least not to close to the crossover glass temperature 
and at a qualitative level, allows us to claim that these extremely long times scales are 
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unrealistic if not too close to the glass transition T g even as far as dynamic fluctuations 
are concerned. 

The ideas discussed in this paper should not only apply to systems that relax 
in a non-equilibrium manner as glasses but also to systems with slow dynamics and a 
separation of time-scales that are kept out of equilibrium with a (weak) external forcing. 
Recently, there has been much interest in the appearance of shear localization, in the 
form of shear bands, in the rheology of complex fluids. Along the lines here described 
it would be very interesting to analyze the fluctuations in the local reparametrizations 
in the fluidized shear band and the 'jammed' glassy band. 

The analytic treatment of mean-field quantum glassy systems follows similar steps 
to the ones presented here. The Schwinger-Keldysh approach replaces the Martin- 
Siggia-Rose one but these formalisms are very similar indeed. The analytic solution to 
the dynamic equations in the limit in which the coupling to the environment is weak 
also uses the fact that the dynamics in the aging regime is very slow. The approximate 
equations then become time-reparametrization invariant. One can then expect that 
similar dynamic fluctuations arise in glassy problems in which quantum fluctuations 
are important. Moreover the proof of global time-reparametrization invariance for spin- 
glasses has been presented directly in the quantum formalism. Novel experimental 
techniques may be apt to study dynamic heterogeneities in glasses when quantum 
fluctuations are important. 

We expect to find similar fluctuations using finite size systems and examining the 
behaviour of the mesoscopic run-to-run fluctuations of the global correlations. These 
may be easier to measure experimentally using mesoscopic systems. 

Importantly enough, global time reparametrization invariance does not develop in 
all models with slow and aging correlation functions. The O(N) ferromagnetic model in 
the large N limit is a case in which global time-reparametrization invariance is reduced 
to just scale invariance [8]. 

Last, but not least, the approach based on reparametrization invariance suggests 
that it may be possible to search for universality in glassiness. A Ginzburg-Landau 
theory for phase transitions captures universal properties that are independent of the 
details of the material. It is symmetry that defines the universality classes. For example, 
one requires rotational invariance of the Ginzburg-Landau action when describing 
ferromagnets. Time reparametrization invariance may be the underlying symmetry that 
must be satisfied by the Ginzburg-Landau action of all glasses. What would determine 
if a system is glassy or not? We are tempted to say the answer is if the symmetry 
is generated or not at long times. Knowing how to describe the universal behavior 
may tell us all the common properties of glasses, but surely it will not allow us to 
make non-universal predictions, such as what is the glass transition temperature for a 
certain material, or whether the material displays glassy behavior at all. This quest for 
universality is a very interesting theoretical scenario that needs to be confronted. 

We have tried to state as clearly as possible the implications of our proposal. The 
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full project is not yet complete since several questions about its limitations remain open. 
A number of issues should be addressed are: 

(i) From a phenomenological point of view, to perform strong tests of the global 
time-reparametrization invariance scenario in molecular dynamic simulations [28~ t |4"4"1 14"5] 
of realistic glassy systems and experiments f!6]-[5T]. More precisely, the triangular 
relations between local coarse-grained correlations and the local fluctuation dissipation 
relations should be analysed to give support or else falsify this conjecture. 

(ii) From an analytic point of view, to derive the effective action for local 
reparametrizations for glassy models with and without quenched disorder. We are 
currently working on this project in collaboration with S. Franz. One idea is to study 
the p spin disordered model with Kac long-range interaction. ANothe one is to study 
the symmetries of the dynamic action associated with the Dean-Kawasaki equation for 
the evolution of the density of a system of particles in interaction. 

(iv) From a mixed analytic and numerical point of view to analyse fluctuations 
in models with global aging dynamics of different type. To this end, one can study 
dynamic fluctuations in simple coarsening systems in finite dimensions. In these cases 
the morphology of domains can be characterized in great detail |43j. We could, in 
principle, understand the fluctuations in the local correlations and linear responses from 
a microscopic point of view. Whether these are similar or different to the ones in other 
glassy problems is still to be established and the outcome of this study could clarify the 
relevance of the value of the effective temperature in determining the characteristics of 
the dynamic fluctuations. 

(v) In the same line as the above, the analysis of fluctuations of elastic lines in the 
presence of impurities could help us understanding the coarsening phenomenon but also 
the role of diffusion that superimposes in these cases to the aging dynamics [T5] . 

These are just a few questions posed by this proposal that are still not answered. 
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